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Abstract 

A pseudoclassical model, reproducing, upon quantization, the 
dynamics of the chiral sectors of the massless spin- 1/2 field theory 
is proposed. The discrete symmetries of the action are studied in 
details. In order to reproduce the positive and negative chiral sectors 
of the particle and antiparticle, we promote the algebra of functions 
on the phase space to a bimodule over the complexified quaternions 
- biquaternions. The quantization is performed by means of strictly 
canonical methods (Dirac brackets formalism) reproducing the Dirac 
equation in the Foldy-Wouthuysen representation in the particle and 
antiparticle sectors and the Weyl equation in the chiral sectors. 
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1 Introduction 



Canonical quantization is the art of establishing a correspondence - following 
Dirac rules - between a classical (or pseudoclassical) model and a quantum 
system. Its universal status is questionable, since one cannot discard the 
possibility quantum mechanics to exist - in some particular cases - on its 
own right, without any relation to some (pseudo-) classical mechanics. On 
the other hand, canonical quantization has proven its efficiency on a vast 
variety of examples. 

The inverse problem - that of finding a (pseudo-) classical system 
reproducing, upon quantization, given quantum dynamics - is also 
interesting. In particular, after this problem being solved, a path integral 
for the quantum system can be obtained. 

The pseudoclassical relativistic particle models [H |2| generally reproduce 
the consistent quantum mechanics [3] of the one-particle and one-antiparticle 
sectors of the corresponding quantum field theory. Two cases have resisted 
the attempts to obtain a pseudoclassical model yelding a consistent quantum 
mechanics upon canonical quantization. These are the massive spin- 1/2 
particle in odd- dimensional space-time (with a given value of spin) and the 
massless particle. In both cases the problem is to obtain a "minimal" model. 
Indeed, a theory, reproducing, for example, the quantum mechanics of a 
spin-1/2 particle with both signs of spin in (2+l)-dimensional space-time, 
can be easily obtained by reduction from the spinning particle model in 
(3+1) dimensions. A common technical problem arises in the both cases: a 
first-class bifermionic constraint which does not admit gauge fixing 

In the present paper we propose a pseudoclassical model for the massless 
spinning particles which is free of the problem. The pseudoclassical action is 
extracted from a path integral in a way similar to that of Refs. jSUHl- in order 
to obtain a canonically quantizable model we had to solve the bifermionic 
constraint problem, as well as the bifermionic constant problem jl]. We 
looked for a model satifying certain condititions. The most important among 
them was to obtain a bifermionic constraint belonging to a set of second- 
class constraints, instead of being a first-class constraint. Also, we looked for 
a model with the simplest possible constraint structure and, in particular, 
leading to the same set of first-class constraints as in the (massless) spinless 
case. 

The canonical quantization of the model reproduces the quantum 
mechanics of the one-fermion state and that of the negative-energy states, 
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corresponding to the antifermion states. 

In order to reproduce the dynamics in the chiral sectors, i.e., the dynamics 
of states of given chirality, we propose a quaternionic generalization of the 
pseudoclassical mechanics. The 'algebra of functions on the the phase space' 
is assumed to be an algebra over biquaternions (a quaternionic bimodule). 
This allows us to develop a pseudoclassical scheme which reproduces, upon 
canonical quantization, the quantum mechanics in the chiral sectors. The 
paper is organized as follows. In Sect. 2 pseudoclassical action is derived 
from a a path integral for the propagator. The discrete symmetries of the 
action are studied in Sect. 3. The Hamiltonian formulation is given in Sect. 4. 
Canonical quantization is performed in Sect. 5. 



2 The path integral 

The propagator S c of the free massless Dirac field satisfies the equation 

ij^S c (x) = -6 4 (x) (1) 

and Feynman assymptotic conditions. In order to get an action, leading to 
a suitable set of constraints, we obtain a path-integral expression for the 
function 

A(x)= 7 5 7°5 c (x). (2) 

It can be expressed, following Schwinger [7j, as a matrix element of an 
operator in a Hilbert space. Let X M , P M be a set of operators, generating an 
irreducible representation of the canonical commutation relations 



P P 

1 ill 1 V 



0. 



X*,P V 



in a Hilbert space, in such a way that 

(x | Pfj, | V) = -id^(x | ip) 

where | x) is an eigenvector for all and ip is an arbitrary vector from a 
suitable dense set in the Hilbert space. Then 



where 



A(iC ut ^in) ( | A 1 3?i n ) 
_ 7 5p0 + 7 5 a i_p. 

P 2 + ie ' 



(3) 
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a 



7°7 (k = 1,2,3), and the weak limit e — > 0+ is understood. We 
omit the term it in the sequel. One can use the Schwinger proper-time 
representation [7] for the inverse operator 1/-P 2 , 



1 

J~2 



dX e 



i\P 2 



(4) 



As to the factor 7 5 P° + 7 5 a • P, it could be represented by an integral over 
a pair of Grassmann variables [Uj, Xi an d X2 , 



7 5 P° - 'fa* Pi = i I dxidx2 exp 



1X1X2 1 7^° + 7^ 



(5) 



Some other possibilities also exist. Let B be some quantity, then 



2tt J c z 2 



(6) 



where C is a curve circling around the origin in the z-plane. Using eqs. Pj). 
© and 

€ij k a j a k = li^fa 1 

in Q, one obtains 



dA i" (x out I exp 
^ Jo Jc z 



i\p 2 



(7) 



+iz( 1 s P + -e ijk P i afa/°)\ |x in >. 



The integrand can be considered as an evolution operator for a quantum- 
mechanical system with Hamiltonian 



H(P, A, z) = -sP 2 - z ( 7 5 P° + -e ijk P i a j a k ) . 



(8) 



A path integral for the matrix element in (|ZJ) is obtained in the standard way 
jHlin]. Using the completness relation 



d 4 x\x)(x\ 
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and introducing 5-functions in s, u and v, one obtains 

(^lexpI-^P,^)]!,,^/^,.../^ 

poo poo 

x / ds\ . . . I dsN S(sn — sjv-i) • • • S(si — so) 
Jo Jo 

x J dui . . . J du N 5{u N — u N _i) . . .5{ui — uq) 



x / dvi 



dv N 5(v N - vn-x) ...5(v x - v Q ) 

(x out | exp{-iH(P, s N ,u N + iv N )A N r} | x N -i) 
(xi | exp{-iH(P,si,ui +w 1 )A 1 t} \ x in ), 



where Y2f=i ^ r = 1- Evaluating the matrix element, 
(x k | exp - iH(P, s k , u k + iv k )A k r \ x k _i) 



I 



(27T)- 



■ cxp 



W 2 * - a*-i) - s fc , M fc + iv k )A k r 



and using integral representations for the 5-functions one obtains 
(x out | exp [ - iH(\, z)] | 

d A XN-l 

-+oo 



= /**.../ 

/■ + OO /■ 

/ dsi . . . / 

Jo Jo 

dui 



(27r) 4 "'y (2tt) 4 



x 



x 



dsN 



+00 



dPsi 
2ir 



dP. 



sN 



X 



X 



00 
+00 



00 
+00 



dh 



N 



dvi 



00 
+00 

00 

oc 



dP, 



hi--- 



00 
00 



du 



N 



/+00 p 
dP ul ... 
-00 J — 



oc 
+00 

■00 
+00 



2tt 



dP, 



hN 



dP 



uN 



00 
+00 



N 

xTexp [pfe- 
fc=i 



dVN 
%k -Kk—1 



dP vl 



OC' 

+00 



dP 



vN 



A k T 



D S k — S fc _i fofe — 

+ -Tsfc 7 r -Oife" 



A k T 



A k T 



+p 



uk' 



A k r 



+ P 



vk ' 



A k r 



-Ti(Pk,s k ,u k + iv k ) A fe r|, 



where u + iv = z, s = A, h = 7 s and T indicates an ordered exponent 
(a-matrices in the factors with greater value of k stay to the left of the ones 
with smaller k). Making the rescaling A = eo/2, Sk = efc/2, P s k = 2P e k and 
taking the limit A^r — > 0, one obtains the path integral 



A c (x out — x- m ) 



4vr Jo 



dec 



dz 

~~2 



C z 



Dx I Dp 



ho 



De I DP P I Dh I DP h I Du / DP, / Dv / DP, 



<'(> 



x T exp < z 



+ P e e + P h h + P u u + P v v 



1 2 

+ -ep 2 + (u + iv) (hp + -eijkpWa*) 



dr 



ft-0=7 5 



The time-ordered exponent in eq. can be expressed as a linear 

combination of antisymmetrized products of a-matrices by means of a path- 
integral technique [21 EH El) with the result 



47T 



Sym 



de r 



dz 

~2 



c z 



%out 



■' t h 



Dx Dp De DP ( 



to 



X 



ho 



Dh / DP h / Du / DP U / Dv / DP, 



x exp < 1 



V(o)+v(i)=e 

P e e + P u ii + P> + P h h + p^ + + 



;i2) 



+ -ep 2 + (u + iv) (hp + 2ie ijk p i ip j ip k ) 



dr - if)(l) ■ V(0) 



0=ol, feo=7 5 



where i/ ? ( r ) is a Grassmann-odd trajectory, # fc are Grassmann-odd 
parameters, and the translationally-invariant measure T>i\) is normalized, 



/ 



T>i\) exp 



i/>(0)+i/>(l)=0 



ift ■ ift dr = 1. 
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Performing the integration over the trajectories in the momentum space p 
one obtains the path integral in the Lagrangian form 



A c (x out ,x in ) = ^Sy m / de 



/ de & — / Dx M (e) / Vij) 

J0 JC z Jx ln J %l>{0)+%l>(l)=G 



eo 



ho 



x De DP e Dh DP h / Du / DP U / Dv / DP V 



UQ 



x exp < % 



■ i]> — —w 2 

2e 



P e e + P u u + P v v + P h h + p^ 
dr-^(l) • V(0) 



0=Ct,/lQ=7 5 



where 



and 



w ° = x° + (u + u> fc = x fc + 2i(u + iv)e klm ijj l ijj m 



M(e) = y Dp exp |^ y ep 2 dr j . 



(13) 



(14) 
(15) 
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3 The action and its symmetries 

The gauge-fixing terms P e e, P u v, P v v, Phh and the symbol-defining term 
—■0(1) • V'(O) apart, the argument of the exponent in (JT3j) is the action. Since 
the complex nature of the variable z = u + iv is not essential, we put v — 
and the action takes the form 

S = I LdT = I ("^ + ^ ' ^) rfr ' ( 16 ) 

where 

w ° = x° + uh, w k = x k + 2iue klm if> l if> m - (17) 

The action ()16|) is reparametrization invariant. It is also invariant under 
space reflection V if the variables are transformed as follows 

V: x k ^-x k , ilj k ^-*jj k } h^-h, u^-u, (18) 

and under the time reflection T', 

T' : x° -> x fc -> ip k ^-ip k , e^e h — >• — /i, it -> it. 

(19) 

The dynamics is invariant under the trajectory- reversal operation 7?., i.e., 




where L is obtained by the following substitutions in L, eq. (JTEJ): 

ft: ^(r) -> x'»{t) = x"(-t), ^*(r) -> ^ fc (r) = -^ fe (-r), 

/i(r) -> fc(r) = -h(-r), u(t) -> u(r) = u(-r), (20) 
e(r) — > e(r) = e(— r) i — ► — i. 

Consequently, the equations of motion are invariant under the 
(pseudoclassical counterpart) of the Wigner time inversion operation 

T = T'K = TIT. 
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4 Pseudoclassical mechanics 



The pseudoclassical mechanics can be formulated in terms of the Z^-graded 
associative algebra A, generated by the canonical variables. For some of 
the non-nilpotent variables we assume the existence of the corresponding 
inverses. The phase space T is a supermanifold 9j. We are not going to 
specify the structure here. However, we refer to the elements of A as to 
functions on the phase space. 

The evolution of a function on the phase space with the evolution 
parameter r is governed by the Hamiltonian equations. The evolution 
parameter r can also be understood as a Grassmann-even generating element 
in a larger algebra A', generated by the canonical variables and r. The 
element r has vanishing Poisson brackets with any function on the phase 
space. 

The constraints are understood as conditions, defining an equivalence 
relation in A. We limit ourselves to the case of r-independent constraints. 
Two elements / and g of A' are considered equivalent if the difference f — g 
is "proportional to the constraints", i.e., if 

f-g = Y,^aG h 

a,i 

where $ a (a = 1 . . . N) are the constraint functions and Fj, are algebra 
elements. Obviously, this define a quotient algebra A (the algebra of 
functions on the constraint suface). Suppose that the set of all constraints 
is a second-class one. If Dirac brackets with respect of this set are used, the 
Hamilton equations define a consistent evolution in the quotient algebra A. 



4.1 Hamiltonian formulation 

The canonical momenta P e , P u , and Ph, conjugate to e, u, and h, 
correspondingly, vanish. The momenta conjugate to tp k are given by 

tt d *- L -ik 
Il k := — — = lip*. 

dtp k 

There are six primary constraints in the theory. The constraint functions are 

$S 1} =P e , ^=Pu, ® { 3 ] = Ph, ^=Iii-l^. (21) 
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and the Hamiltonian if ^ is given by 

# (1) = ~\ep 2 - u (h Po + 2ie klm p k ^ m ) + A a $«, (22) 

where X a are Lagrange multipliers associated with the constraints. The 
constraints $ 4 j can be eliminated. Under a change of the fermionic variables, 

1> k = \ + V k ), n fc = i - V k ) (23) 
one finds that the new variables have the following Poisson brackets: 

{CC 1 } 2M, r {T H ,T h } = -2i5 ij , {&,!&} = 0. (24) 

and $4^ = £\ The conservation of the constraint functions implies 
A 4 * = and the terms in linear in £ fc , vanish. On the other hand, the 
Poissonn brackets with all the variables of the terms in bilinear in £ fc , 

vanish on the constraint surface and can be omitted. The Hamiltonian takes 
the form 

# (2) = ~\e P 2 - u (hp + % -e iik Hrf] + E ^ '■ ( 25 ) 

V / a =l 

The conservation of the primary constraints an d yield three 

secondary constraints, 

$S 2) =p 2 , $< 2) = + \ ®f ] =u. (26) 

The consistency conditions on the constraint functions $ 2 , $3 imply 
A 3 = A 2 = 0. The only Lagrange multiplier remaining unfixed is A . 
One can separate the constraints into a set of first-class constraints, 

^ ] =P e , $S 2) =p 2 , (27) 

and a second-class set of constraints, 

$ 2 1} = P u , $f = u, (28) 

= P h , Of® = hp + l -e ijk pWv k - (29) 
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The set of the first class constraints is exactly the one appearing in the 
theory of scalar massless particles. Only one constraint function depends 
on the Grassmann-odd variables rj k and it belong to the set of second-class 
constraints. 

In this stage $2 i s already contained in a set of second-class constraints, 
that is, without using any gauge condition we get to solve the bifermionic 
constraint problem [Ij. 



4.2 Gauge fixing 

With the aim to use the physical time x° as evolution parameter, we impose 
the gauge condition 

$f = 0, where $f = a; - r. (30) 

The constraint function $f depends on r. In order to obtain a r-indepent 
set of constraints we make a canonical transformation ^Uj- The generating 
function is 

W = x% - rp' , 

so that the relation between the new variables x ,fl , p'^ and the old ones, x ,tJl , 
P'u,, is given by 

X '0 = x 0- T) x ' k = x \ p' ii =p ii . (31) 
The constraint function $f takes the form 

$f = x'° (32) 
and the new Hamiltonian is given by 

# (3) = -Po - \ep 2 - u (hp + ^UjkP'rfrn + A 1 ^. (33) 

The conservation condition for the constraint $f implies 

$2^ = where $^ = ep + uh + 1. (34) 

The consistency condition for this constraint yields A 1 = 0. One can replace 
the set of constraints obtained by the equivalent set <3>, 

$1 = P e , $ 2 = ep + 1, $ 3 = P u , $ 4 = u, 

$ 5 = p 2 , $ 6 = x'°, $ 7 = P h , $ 8 = hp + y ijk p i V j V k - (35) 
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This is a set of second-class constraints as a direct calculation proves. 

There are seventeen fundamental variables at this stage and eight 
constraint equations. Apparently, nine variables should be considered 
independent and the rest of the varibles can be expressed in terms of the 
independent ones. However, one of our constraints is quadratic. 



4.3 Dirac brackets 

The nonzero Dirac brackets of the canonical variables with respect to the set 
of constraints $, eq. ()35|). are 

f j ~i Pi f i ~) Pk 

{^W) = ^*pV- (36) 
v ' Po 

Only the term — po in the Hamiltonian does not vanish on the constraint 
surface. Then the Hamiltonian 



#phys = ~P0 (37) 

can be used instead oh . The evolution of the quotient algebra elements 
is consitently defined by the Hamiltonian equations. However, the Dirac 
brackets do not possess canonical form. This was to be exepcted, because 
the variables p^ are not independent: they are related by the constraint 
equation $ 5 = 0. 



4.4 Solutions in the complex algebra 

If the pseudoclassical algebra is just a C-module, only two solutions exist for 
the equation <3> 5 = 0, namely, 

Po = -(Vtf, (38) 

where ( = ±1 is the energy sign. Putting ( = 1 one obtains, upon 
quantization, the quantum mechanics of the spinning massless particles with 
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both signs of helicity. The second solution yields the quantum mechanics of 
the negative-energy states with both signs of helicity The chirality h is equal 
to helicity, 

h = —L=(e ijk pWr] k . (39) 
2 v P 

and is a conserved quantity. Indeed, its Dirac bracket with the Hamiltonian 

k 

vanishes. 



H ( = CVp 2 



4.5 Quaternionic pseudoclassical mechanics 

Since our aim is to reproduce, by strictly canonical methods, the quantum 
mechanics of the chiral sectors, we assume that the pseudoclassical algebra 
A is a (B, B)-bimodule, where B is the algebra of the biquaternions. Further, 
we postulate the following relation between the left and right multiplication 
of the algebra elements by quaternions HI. Let {1, J 1 , J 2 , J 3 } be a basis in 
H, such that 

J*jS = -5V + eijk J k . (40) 

Then 

J7] \ v k J* if i = k 1 > 

The rest of the generating elements of A, as well as the evolution parameter 
r, commute with J k . These rules are invariant under a change of the basis 
in H 

{1, J\ J 2 , J 3 } - {1, J' 1 , J' 2 , J' 3 }, 

if J' 1 , J' 2 , J' 3 satisfy eq. fljOj) , because the quotient algebras arising when the 
two bases are used, are isomorphic. 

Equation p 2 = possesses new solution in the extended algebra, in 
addition to the solutions (|3*Sjl. These are 

where J fc G H satisfy (J4(J|) . The constant x — =tl is the (pseudoclassical 
counterpart) of helicity (multiplied by two). Then 

Po 1 = ~\iJ k Pk- 
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Because of the relation postulated flU)) between the left and the right 
multiplication by quaternions, the constraint equation <3> 8 = is ambiguous 
without specifying the order of the factors in its right-hand side. Assuming 
that h stays to the left of po, one obtains 

The quantity h is conserved, since its Dirac bracket with the Hamiltonian 

H = i X J k p k (42) 

vanishes. 

5 Quantization 

According to Dirac rules, the operators x k , pk, and f) h , corresponding to the 
independent variables x k , p k , and r] k , must satisfy the (anti-) commutation 
relations 

[x\pj]_ = i5), (43) 
[f,\rf] = 25*. (44) 



5.1 Positive- and negative-energy sectors 

Let TC xp be a Hilbert space in which an irreducible representation, generated 
by x\ pi, of the Heisenberg algebra is given. Let H v be a unitary space in 
which an irreducible representation, generated by rf , of the Clifford algebra 
is given. Then in H = Ti. v <8> 7~t xp an irreducible representation of the Lie 
superalgebra is given. If the solution (jHEJ) is chosen, then 

Po = -cV¥- 

The energy sign operator 




is a number, and, therefore, it commutes with all the canonical operators. 
The rest of the dependent operators, h and e, are expressed in terms of the 
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independent ones, 



h = ^h e vkP if l j v\ e = (45) 
2a/p vP 



according to eqs. (|35p. and the Hamiltonian is given by 

H c = (V¥- (46) 

The Schrodinger equation reads 

" () 1 ■(x°)) = CV¥\<P(x°)), (47) 



where | 4>{x )) is a time-dependent vector in 7i. An explicit realization is 
given in the Hilbert space of the (two-component, square-integrable) spinor- 
valued functions 0(x), where one can put 

d 

x k = x\ p k = -i—, v k = a k , (48) 

o k being the Pauli matrices. Taking the direct sum of two copies of this 
Hilbert space and putting ( = 1 for the first copy and ( = — 1 for the second, 
one obtains 

C = P = 7°, V k = 7 5 7V- 

The Schrodinger equation for this system is the Dirac equation in the Foldy- 
Wouthuisen form, 

i^{x)=(3^^{x). (49) 
and the chirality operator h reads 

h = ITJ^^^ffrf = ^27 5 7 fc Pfc- 

2 VP 2 VP 

It anticommutes with the operators fj k (and commutes with x h , pk), while 
the energy-sign operator ( commutes with all the canonical operators. 
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5.2 Chiral sectors 

Let a = a + di J % G B be, where do, aj G C. In biquaternionic algebra we can 
define three conjugations, namely 

a =a + a i J, 

CL = CLq — QiJ Z , 

a ' = a = a — J . 

These are complex, quaternionic and hermitean conjugations, respectively. 
With these operations, we can define biquaternionic scalar product as 

(a, b) = Re{a j b). 

Following [TT] we have 

Re(z) = ^(z + z), 

and also (zw) = wz, w — w. Then we can rewrite the scalar product as 

(a,b) = ^(a*b + ba*). (50) 

The operators r] k , acting in B and defined by 

rj k = iJ k , 

satisfy the anticommutation relations (|44|) and posses the same commutation 
properties with J k as the pseudoclassical variables r/ k do. The representation 
of the Clifford algebra in B, generated by these operators is reducible and 
decomposes into two irreducible representations. The subspace B + , spanned 
on the vectors 

|l) = -L(l + zJ 3 ), |2) = i=(j 2 + ^). 

is one of the invariant subspaces. The representation of the Clifford algebra 
in B + , generated by f) k is irreducible. Let an irreducible representation, 
generated by x k and pk, of the Heisenberg algebra be given in some complex 
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Hilbert space 7i xp . An irreducible representation of the canonical (anti-) 
commutation relations is defined in the tensor product TC = B + <g) 7i xp in a 
natural way. The Hamiltonian is given by H = ixJ k pk, according to eq. (|42|) . 
Let | x) be the basis of the coordinate eigenvectors in TC xp , 

d 

x fc | x) =x fc | x), (x I p k I 0) = -i— (x I 0), (x| x) = 5 3 (x-x'). 

A basis in 7i is given by the vectors | a, x) =| a)® x) (a = 1,2). A 
state vector | ip) in Tt is represented in this basis by a two-component wave 
function ip, 

(l,x| V) " 
(2,x| V) 



^(x 
Since 



(l,x | r/ fc | V) \ _ k 



( 2 ,x|^i^) ;- a '^ x )' 

the time dependent wave function satisfies the Weyl equation 
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